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Spontaneous magnetization of quark liquid is examined. It is pointed out that 
quark liquid has potential to be ferromagnetic at rather low densities. 

1 Introduction 

Pulsars are rotating neutron stars emitting radio waves, X-rays or gamma 
rays. Ordinary radio pulsars have a magnetic field of O(10 12 ~ 13 )G, which 
causes various radiations. The origin of such strong magnetic field is still an 
open problem. Recently a new type of neutron stars, called magnetars, has 
been proposed to explain the observational data_pn pulsars, which should have 
an extraordinary magnetic field of O(10 15 ) G El. There are reported several 
magnetar candidates so far for anomalous X-ray pulsars (AXP) and pulsars 
associated with soft-gamma-ray repeaters (SGR). 

There has been a naive working hypothesis to understand the magnetic 
field in neutron stars; if the magnetic flux of a main sequence star is conserved 
during its evolution, the decrease in radius leads to an increase in the magnetic 
field. For example, the sun, a typical main sequence star, has a magnetic field 
of O(10 3 )G with the radius R ~ 10 10_11 cm. By squeezing the radius to 10 6 cm 
for neutron stars we have O(10 11_13 )G, which is consistent with observations 
for radio pulsars. However, if this argument is extrapolated to explain the 
intense of the magnetic field for magnetars, their radius should be O(10 4 )cm, 
which is much less than the Schwartzschild radius of neutron stars with the 
canonical mass M = 1AM Q , R Sc h = IGMjc 2 = 4x 10 5 cm. 

These observations seems to enforce our reconsideration of the origin of the 
magnetic field in neutron stars. Since there is a bulk hadronic matter beyond 
the nuclear density (rig ~ 0.16fru ) inside neutron stars, it should be inter- 
esting to consider the hadronic origin of the magnetic field; ferromagnetism 
or spin-polarization of hadronic matter may give such magnetic field. Unfor- 
tunately there has been little suggestion about the possibility of spontaneous 
magnetization of hadronic matter. We consider here the possibility of fer- 
romagnetism of quark liquid interacting with the one-gluon-exchange (OGE) 
interaction^. 

One believes that there are deconfinement transition and chiral symmetry 



restoration at several times the nuclear density, while their critical densities 
have not been fixed yet. One interesting suggestion is that three-flavor symmet- 
ric quark matter (strange quark matter) ground or above the nuclear density 
may be the true ground state of matter 0. If this is the case, strange quark 
stars, where quarks occupy almost whole the inner region of stars, can exist 
in a different branch from the neutron-star branch in the mass-radius plane. 
Otherwise quark matter may exist in the small core region of neutron stars. 
We shall see our results should give an origin of the strong magnetic field in 
the context of strange quark-star scenario. 



2 Ferromagnetism of quark liquid 

Quark liquid should be totally color singlet (neutral), which means that only 
the exchange interaction between quarks is relevant there. This may remind 
us of electron system with the Coulomb interaction in a neutralizing positive 
charge background. In 1.929 Bloch first suggested a possibility of ferromag- 
netism of electron systema He has shown that there is a trade off between the 
kinetic and the exchange energies as a function of density, the latter of which 
favors the spin alignment due to the Pauli principle. This was a beginning of of 
the concept of itinerant magnetism. In the following we discuss the possibility 
of ferromagnetism of quark liquid on the analogy with electron gas. 

It is to be noted that there is one big difference between them; quarks 
should be treated in a relativistic way. The concept of the direction of spins is 
not well defined in relativistic theories, while each quark has two polarization 
degrees of freedom. Here we define the spin-up and -down states in the rest 
frame of each quark. Then the projector onto states of definite polarization is 
given by 

P{a)= l -{1 + 1 ^) (1) 

with the 4-peudovector a, 

— k(C ' k) a°=<^ (2) 



m q (E k +m q )' 

for a quark moving with the momentum k = (E k , k) i. The 4-peudovector a 
is reduced into the axial vector £ (|£| = 1) in the rest frame, which is twice 
the mean spin vector in the rest frame. Actually if we choose C along the z 
axis, £ = (0, 0, ±1), we can see each value corresponds to the spin-up or -down 
state. The mean value of the spin is given by 

1 m q ( k(C • k) 



2E k \ C+ m q {E k +m q )) ' 
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Finally the projection operator P(a) gives the polarization density matrix p, 



p(k,Q = + m q )P(a), P(a) = ±(1 + 75 #. (4) 

The exchange interaction between two quarks with momenta k and q is 
given by 

/kc.qC = -j^^-M^,^,. (5) 

•MkC.qC' is the usual Lorentz invariant matrix element, and is evaluated with 
the help of the polarization density matrix (Q) 

MkC.qC = 5 2 ^tr(A a /2A Q /2)ltr [ 7M p(fc, ChM?, C)] (6) 

2 1 

= g 2 -[2m 2 q -k-q-m 2 a-b} - _ , (7) 

where the 4-pseudovector 6 is given by the same form as in Eq. (|^) for the 
momentum q. 

Although the vector £ of each quark may point in a different direction on 
the two dimentional sphere S 2 , we assume here it along the same direction, say 
z axis. The exchange energy for quark liquid is then given by the integration of 
the interaction (^|) over the two Fermi seas with the spin-up and -down states; 
eventually, it consists of two contributions, 

_ non-flip , flip /o\ 
tea; — t ex T t ex ■ {O ) 

The first one arises from the interaction between quarks with the same po- 
larization, while the second one with the opposite polarization. The non-flip 
contribution is the similar one as in electron gas, while the flip contribution 
is a genuine relativistic effect and never appears in electron gas. We shall see 
that this relativistic effect leads to a novel mechanism of ferromagnetism of 
quark liquid. 



3 Examples 

We show some results about the total energy of quark liquid, e t ot — tkin + £ex 5 
by adding the kinetic term e^m- Since gluons have not the flavor quantum 
numbers, we can consider one flavor quark matter without loss of generality. 
Then quark number density directly corresponds to baryon number density, if 
we assume the three flavor symmetric quark matter as mentioned in §1. 
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Figure 1: Total energy of quark liquid as a function of the polarization parameter for densities 
n q = 0.1,0.2fm- 3 . 

There are two parameters in our theory: the quark mass m q and the quark- 
gluon coupling constant a c . These values are not well determined so far. In 
particular, the value of quark mass involves subtle issues; it depends on the 
current or constituent quark picture and may be also related to the existence 
of chiral phase transition. Here we allow some range for these parameters and 
take, for example, a fiducial set, m q = 300MeV for strange quark and a c = 2.2, 
given by the MIT bag modelQ. In Fig. 1 two results are presented as functions 
of the polarization parameter p defined by the difference of the number of the 
spin-up and -down quarks, n+ — n~ = pn q . The results clearly show that 
the ground state should be ferromagnetic for lower density, while it is in the 
paramagnetic phase for higher density. The phase transition is of first order 
and its critical density is around n q ~ 0.16fm -3 in this case, which corresponds 
to the nuclear density for flavor symmetric quark matter. Note that there is 
a metastable ferromagnetic state (the local minimum) even above the critical 
density. This ferromagnetic phase is a spontaneously symmetry broken state 
with respect to the rotational symmetry: the order parameter is the mean 
value of C, (C)j an d symmetry is broken from G = 0(3) to H = 0(2) once (£) 
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Figure 2: Density dependence of three quantities. The negative regions for Se and r) indicate 
the ferromagnetic phase and the existence of the metastable ferromagnetic state. 

takes a special direction on S 2 . 

Magnetic properties of quark liquid are characterized by three quantities, 
Se, x and 77; 8e = e to t{p — 1) — e tot(p = 0), which is the measure for ferromag- 
nctism to appear in the ground state. For small p« 1, 

etot - e tot (p = 0) = X ^P 2 + 0(p 4 ). (9) 

X is proportional to the magnetic susceptibility and plays an important role 
if the phase transition is of second order. In our case it is less relevant since 
the phase transition is of first order. Finally, i] = de to t/dp \ p =i, which is 
the measure for metastability to to exist. In Fig. 2 the density dependence of 
three quantities arc given for a fiducial set of parameters. We can see that 
ferromagnetic phase is the ground state below rig, while the metastable state 
is possible up to rather high densities. 

Finally we present a phase diagram in the m q —a c plane for n q = 0.3fm~ 3 , 
which corresponds to about twice the nuclear density for flavor symmetric 
quark matter. The region above the solid line shows the ferromagnetic phase 
and that above the dashed line indicates the existence of the matastable state. 
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Figure 3: Phase diagram in the mass (m q )- the coupling constant (a c ) plane. <5e in the 
nonrelativistic calculation is depicted for comparison (the dotted line) . 

For massive quarks with the large mass, which may correspond to the current 
s quarks or the constituent quarks before chiral symmetry restoration, the 
ferromagnetic state is favored for small coupling constant due to the same 
mechanism as in electron gas. The ferromagnetic state is favored again for 
light quarks with small mass, which may correspond to the current u, d quarks, 
while the nonrelativistic calculation does not show such tendency. Hence this 
is due to a genuine relativistic effect, where the spin- flip interaction plays an 
essential role. 

4 Summary and Concluding remarks 

We have seen that the ferromagnetic phase is realized at low densities and 
the metastable state is plausible up to rather high densities for a reasonable 
range of the QCD parameters. If a ferromagnetic quark liquid exists stably or 
metastably around or above nuclear density, it has some implications on the 
properties of strange quark stars and strange quark nuggets. They should be 
magnetized in a macroscopic scale. For quark stars with the quark core of r q , 
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simply assuming the dipolar magnetic field, we can estimate its strength at the 
surface R, 




amounts to order of 0(1O 15 ~ 17 )G for r q ~ 0(R) and n q = O(0.1)fm~ 3 , 
which should be large enough for magnetars, using the quark magnetic moment 
Hq ~ Hn([J-n muclear magneton~ 5 x 10~ 24 erg • gauss -1 ) for massive quarks 
and 10 2 /Uat for light quarks. Hence it might be interesting to model SGR or 
AXP using our idea. 

We have found that ferromagnetic instability is feasible not only in the 
massive quark system but also in the light quark system: the spin-nonflip 
contribution is dominant in the nonrelativistic case as in electron gas, while a 
novel mechanism appears as a result of the large spin-flip contribution in the 
relativistic case. 

Our calculation is basically a perturbative one and the Fermi sea remains 
in a spherical shape. However, if we get more insight about the ferromagnetic 
phase, we must solve the Hartree-Fock equation and thereby derive a self- 
consistent mean-field for quark liquid. Moreover, we need to examine the long 
range correlation among quarks by looking into the ring diagrams, which has 
been known to be important in the calculation of the susceptibility of electron 
gas. 
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